clc;
syms x y;
f=(x^2-2*x)*exp(-x^2-y^2-x*y);
H=simplify(hessian(f,[x,y]))
H1=simplify(hessian(f,[x,y])/exp(-x^2-y^2-x*y))